Giant acceleration in slow-fast space-periodic Hamiltonian systems 
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Motion of an ensemble of particles in a space-periodic potential well with a weak wave-like per- 
turbation imposed is considered. We found that slow oscillations of wavenumber of the perturbation 
lead to occurrence of directed particle current. This current is amplifying with time due to giant 
acceleration of some particles, ft is shown that giant acceleration is linked with the existence of 
resonant channels in phase space. 
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In the recent years considerable interest was devoted to 
the ratchet effect, generation of directed particle current 
in the absence of any biased forces. This phenomenon 
is relevant for wide range of applications, including con- 
trolled photocurrents in semiconductors [H, 0| , motion of 
cold atoms in optical lattices 0, 0, [E| , transport of passive 
tracers in meandering jet flows in the ocean [6|], biolog- 
ical and chemical systems (see Q for a comprehensive 
review). 

The ratchet effect in space-periodic Hamiltonian sys- 
tems is associated with the asymmetry of chaotic region 
in phase space [1, |T(| ■ Recently it was reported about 
new type of Hamiltonian ratchets, in which a periodic 
potential is subjected to a sum of external forces which 
are periodic in time and in coordinate [Til ] . Each of the 
forces induces strong but local chaotic diffusion in certain 
areas of phase space. This effect is achieved by means of 
resonance between temporal and spatial oscillations of 
the perturbation imposed. This resonance is asymmetric 
in momentum space, that provides asymmetry of cross- 
ing the separatrix and occurrence of directed transport. 
Combining such forces, we can make finite motion un- 
stable for all range of the particle energy. By this way, 
even a weak perturbation of the potential can activate 
ballistic current of particles with the lowest initial ener- 
gies. A similar effect was used in [l2[ in order to produce 
surfatron acceleration. 

In the present Letter we demonstrate the mechanism 
providing simultaneously generation and giant acceler- 
ation of directed current by means of a weak external 
perturbation. Possibility of giant acceleration arises due 
to adiabatic variations of the perturbation. 

Consider an ensemble of non-interacting unit-mass 
point particles driven by a wave-like external force. Mo- 
tion of each particle is described by the Hamiltonian 

P 2 

H = — — cos x + e cos(kx + vt), (1) 

where e <C 1, and wavenumber of the perturbation k is a 
slowly- varying parameter 

jfe = fe(l + acosfif), \a\ < 1, 0<1. (2) 



Physically this condition corresponds to slow libration 
of the external force with respect to axis x. Particle 
trajectories obey the Hamiltonian equations 

x = p, p = — sin x + ek sin 0, (3) 

where we denote the perturbation phase kx + vt as <f>. 
Hereafter we shall consider the case when the parameters 
k and v have sufficiently large positive values, so that 
the inequality d<fi/dt ^ 1 is satisfied along a particle 
trajectory, except for some small resonant region, where 

— = kp - kaQx sin fit + v ~ 0. (4) 
at 

Outside this region particle dynamics is close to inte- 
grate and can be described using the averaging tech- 
nique [l3j]. According to (HJ, the resonant area in phase 
space is located along the line given by the following 
equation 

v aQ sin \P 

Pres = ~ k(l + acos*) + l + acosv^' { ) 

where VP = fit. In order to describe the motion inside 
the resonant region we derive, using ([3]) and (|4]), the 
"pendulum-like" equation for (f> 

4>- ek 2 sm(f> + f{x,p, *) = 0, (6) 

where f(x,p, 'J) is treated as a slowly-varying parameter 
given by the following equation 

f(x,p,t) = akQ(2psin^ + fix cos*) + fcsinx. (7) 

Equation © corresponds to the Hamiltonian of the fol- 
lowing form: 

H{4>, 0) = l -4> 2 + <l,f(x, p, *) + e~k 2 cos 0. (8) 
If the inequality 

\f(x,p,*)\<sk 2 (9) 

is satisfied, then the phase portrait corresponding to the 
Hamiltonian ([5]) contains an oscillating resonant region 
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in 4>-<p space, bounded by a separatrix. Falling into reso- 
nance (0]) corresponds to crossing the separatrix and en- 
tering this region. Probability of entering depends on the 
area of the resonant region and increases with increasing 
the difference between the left- and right-hand sides of 
(J9j) . A particle spends any time inside the resonant zone 
and then leaves it with strongly increased or decreased 
energy. Following [Til . [TH [l6| , we derive an approximate 
formula for the energy jump 
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ek 2 cos < 



(10) 



where < <f> < 2tt, /*, p* and cf>* are values of /, p 
and 0, respectively, when the trajectory hits the resonant 
region. AE depends extremely on </>*, therefore multiple 
recurrences to the resonant area cause chaotic diffusion 
in phase space [H EE El- 
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FIG. 1: Function g(x, *) with x = -1000, x 
x = -50000, and x = -100000. 
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With fixed values of k and v 1 the resonant condition 
(|3|) has the simplest form 



kp + v ~ 0. 



(11) 



It should be emphasized, that inequality ((9]) is fulfilled 
only if a particle is not far from an extremum of the 
unperturbed potential. Hence the condition (jlip should 
be replaced by the following one 18]: 



p(E res , x = itl) = Pr. 
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(12) 



where I is integer. Using (fT2")) we find the resonant values 
of energy 



(13) 



Equation (TT5|) determines locations of the chaotic layers 
in energy space [HI [3 ■ If the chaotic layer induced 
by resonance (fT2")) coalesces with the near-separatrix 




4000 6000 

t (a. u.) 



8000 10000 



FIG. 2: (a) Mean coordinate, (b) mean momentum and (c) 
momentum variance as functions of time. 



chaotic layer, then the chaotic sea formed has much larger 
width in the lower half-plane of phase space than in the 
upper one [ll|. It follows from the asymmetry of the 
condition (|12[) in momentum space and implies the preva- 
lence of chaos-induced particle flights towards x = -co. 

It is natural to suggest that adiabatic variation of the 
resonant momentum (|12p leads to a gradual displacing 
of areas of instability in phase space. If the timescale of 
diffusive mixing inside the chaotic areas is much smaller 
than the timescale of changing the resonant value of mo- 
mentum (1121) , then these areas play the role of dynamical 
trap s for particles, so-called stochastic layer traps (SLT) 
[20l |2~0 ] . Consequently displacement of a chaotic layer in 



3 



s 

a 

a. -10 

-15 



a) 



-20 

-4000 -3000 
10 




-2000 -1000 

x (a. u.) 



1000 



3 

a 



10 



-20 



-30 




-4000 -3000 -2000 -1000 

x (a. u.) 



1000 2000 



FIG. 3: Particle distribution in phase space at (a) t = 1300, 
(b) t = 1400. The resonant channel is marked by the line. 



energy space can be followed by increasing or decreasing 
of mean energy of particles belonging to it. As it will be 
shown in this Letter, a rather complicated situation oc- 
curs if wavenumber of the perturbation varies according 
to the law ©. 

Substituting (0 into 0, we obtain the expression for 
the criterion (|9|) on the resonant line 



afi 2 (-acos 2 * + cos* + 2a) 



(1 + a cos*) 2 
2aVlv sin * 



SIM 



k(l + a cos*) 2 



(14) 



< ek. 



This inequality holds if sin x ~ and, subsequently, x ~ 
irl, where I is integer. When skipped the term ~ sinx, 
we can rewrite the criterion (|14p as follows: 



j(x, *) 



aft 2 (-a cos 2 * + cos * + 2a) 



(1 + a cos*) 2 
2aSl^sin* 



k(l + a cos*) 2 



(15) 



ek < 0. 



Figure [T] represents the function g(x, *) with different 
fixed values of x. According to this figure, the criterion 
(P3J) is satisfied with \x\ < 50000 for the large intervals 
of *, centered at 27rm, where m = 0,1,2,.... This im- 
plies existence of those trajectories which, being passed 



the resonance area at once, will visit the resonant area 
repeatedly on the subsequent cycles of pendulum, till the 
slowly- varying phase * remains close to 2irm. Such par- 
ticles move along the lines described by $5§, towards the 
point x — when sin* < and from it when sin* > 0. 
The latter ones are capable to perform ballistic flights 
with increasing velocity. 
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FIG. 4: The same as in Fig. [2] at (a) t = 3200, (b) t = 5200, 
(c) t = 9200. 



Occurrence of such ballistic flights is confirmed by nu- 
merical simulation. We computed evolution of the ensem- 
ble of particles, initially distributed with gaussian prob- 
ability density 



p(x, p,t = 0) 



1 



2tT<J 0x (70p 



exp 



P 



(16) 
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where er 0x = (r 0p = 0.1. The parameters of the perturba- 
tion we used are the following: e = 0.04, k Q — 12, v = 4, 
a = 0.75, £1 = 27r/1000. Figure [2] represents the tem- 
poral dependence of mean coordinate, mean momentum 
and variance of momentum. It is shown that there oc- 
curs a particle flux directed towards t — > — oo. The mean 
momentum grows nonmonotonically and abrupt acceler- 
ations are alternating with abrupt slowings down. Ac- 
celeration takes place when the slow phase 'J is close to 
27TTO, that agrees with our analysis. Each act of accel- 
eration is followed by step-like increasing of momentum 
variance. It should be emphasized that momentum vari- 
ance is much larger than mean momentum, that indicates 
the presence of particles with very high velocities. Fig- 
ure [3] shows that accelerating particles form jets along 
the resonant line ((SJl, which is cut according to the cri- 
terion (|15p. The first significant jet becomes apparent 
at t ~ 1300. Accelerating particles follow the resonant 
line until t ~ 1400 and then leave the resonant zone. It 
should be noted that strongly accelerated particles never 
return into the resonant zone again. Forming of later jets 
is demonstrated in Fig. [H where instantaneous particle 
distributions at t = 3200, t = 5200 and t = 9200 are pre- 
sented. Evolution of the particle cloud is also presented 
in the media files, which are available at [22 |. 

We can distinguish three stages of evolution of the 
particle ensemble. At the first stage directed current 
is activating. Until a particle is not far from x = 0, 
the location of the resonant zone is determined by the 
formula (|lip. i.e. by the instantaneous value of the ra- 
tio v/k(t). It infers that the initial particle cloud is 
placed inside the chaotic layer caused by resonance with 
Pics = —vjk(t = 0) = —4/21. Adiabatic decreasing of k 
displaces this layer to the separatrix. The time of diffu- 
sive mixing inside the chaotic layer is much smaller than 



27r/f2, so that the particle cloud follows the chaotic layer. 
At $ — ir chaotic layer caused by resonance (0| merges 
into the near-separatrix chaotic layer, that leads to oc- 
currence of ballistic particle current towards x = — oo. 

The second stage starts when the particle cloud be- 
comes enough wide and some particles are capable to fall 
into the resonant channels described by $5§. This stage is 
characterized by fast growth of momentum variance due 
to events of giant acceleration. Note that some particles 
turn around and then perform ballistic flights in the op- 
posite direction. Nevertheless, number of particles accel- 
erating in the direction x — > — oo is much larger, therefore 
the turned particles give negligible contribution into the 
resulting particle flux. Since the resonant channels have 
finite length, one can call the zone where they exist as 
the accelerating zone. 

The third stage is not presented in the figures. It 
starts when the particle cloud becomes very wide and 
only negligible fraction of a particle ensemble remains 
within the accelerating zone. At this stage momentum 
variance achieves saturation and stops increasing. 

In conclusion, in this work we demonstrated the effect 
of giant particle acceleration in the simple space-periodic 
Hamiltonian system subjected to a slowly-modulated ex- 
ternal force. The effect arises from the specific topology 
of resonance ^ in phase space, which permits capturing 
of a particle into the accelerating channel. 
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